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Abstract. Evidence is presented that the WIND spacecraft observed particle and 
field signatures on October IS- 19, 1995 due to reconnection near the footpoints of 
a magnetic cloud (i.e., between 1 and 5 solar radii). These signatures include (i) 
an internal shock traveling approximately along the axis of the magnetic cloud, (ii) 
a simple compression of the magnetic field consistent with the footpoint magnetic 
fields being thrust outwards at speeds much greater than the solar wind speed, (iii) 
an electron heat flux dropout occurring within minutes of the shock indicating a 
topological change resulting from disconnection from the solar surface, (iv) a very 
cold 5 keV proton beam and (v) an associated monochromatic wave. We expect 
that, given observations of enough magnetic clouds, Wind and other spacecraft will 
see signatures similar to the ones reported here indicating reconnection. However, 
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these observations require the spacecraft to be fortuitously positioned to observe 
the passing shock and other signatures and will therefore be associated with only 
a small fraction of magnetic clouds. Consistent with this, a few magnetic clouds 
observed by Wind have been found to possess internal shock waves. 

1. Introduction 

Magnetic clouds, a subset of interplanetary ejecta characterized by strong mag- 
netic fields which exhibit a smooth, large rotation and low proton temperatures, are 
currently of great interest in part because of their association with coronal mass 
ejections (CME’s) and because of their magnetospheric impact [Burlaga et al ., 1981; 
Laakso et al., 1998; Moore et al., 1999]. Magnetic clouds typically expand at about 
half the Alfven speed and, due to their adiabatically decreasing magnetic field, are 
expected to disappear somewhere between 2 and 12 AU [Osherovich et al., 1993]. 
With the fleet of satellites at 1 AU resulting from the International Solar-Terrestrial 
Physics (ISTP) program, magnetic cloud observations from solar origin through ge- 
omagnetic effect are now possible [e.g., Fox et al., 1998]. 

The large-scale geometry of magnetic clouds is well- described by a force-free 
magnetic field represented by a set of helical field lines confined to a flux tube 
(flux rope) [ Burlaga , 1995; Wu et al., 1995; Low and Hundhausen, 1995] or flux 
tubes (flux ropes) [Osherovich et al., 1999]. However, there appear to be smaller- 
scale physical processes occurring within magnetic clouds [Christon et al., 1998; 
Takeuchi et al., 199S]. 

Because magnetic cloud footpoints are frequently still attached to the Sun, we 
might expect to observe signatures on magnetic cloud field lines of time-dependent 
processes in the lower corona. One such time-dependent process is magnetic re- 
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connection which effects the disconnection of the footpoints [ Gosling et al., 1995]. 
Indeed, Gosling [1990] posits this as the basic process by which interplanetary flux 
ropes are formed. 

2. The October 1995 Cloud 

One example of a well-studied magnetic cloud was observed by the Wind space- 
craft on October 18-19, 1995. An overview of the magnetic field and plasma pa- 
rameters during this time is shown in Figure 1. The first three panels display Wind 
magnetic field data from the Magnetic Field Investigation (MFI) [ Lepping et al ., 
1995] with the upper panel showing the magnetic field magnitude in nanoteslas on 
a logarithmic scale (see Fig. 1 of Lepping et al. [1995] for a linear scale). The next 
two panels show in degrees the magnetic field angle out of the ecliptic plane, 6, and 
in the ecliptic plane, (j> (with 0° pointing towards the Sun). The interior of the cloud 
is indicated by the two headed arrow between the vertical lines. Although there is 
some ambiguity about the location of the back of the cloud, it does not bear on this 
analysis, so we have only indicated one of the possible boundaries [ Lepping et al., 
1997], 

This cloud has been fit using a force-free “constant a" helical flux rope model 
resulting in cylindrical Bessel functions: 

B^ = B 0 Ji(ar) (1) 

B z = BoJo(ar), (2) 

where B $ and B : are the azimuthal and axial fields, respectively, Bo is the central 
field, r is the distance from the cloud axis, and a is a constant relating the current 
and magnetic field (qB = V x B) [Priest, 1987]. 
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The fitting result was awarded a rating of “1: excellent” on the subjective 
Lepping quality scale. This scale includes the ranks 1 for "'excellent”, 2 for “good”, 
3 for “poor” and “cl” for “cloud-like” indicating that the helical flux rope fitting is 
unlikely to be successful for any reasonable model. Of 34 magnetic clouds identified 
by Lepping using Wind data from 1995-1998 (including the October 1995 cloud), 
13 were rated “excellent” (1), 15 were rated “good” (2), 5 were rated “poor” (3) 
and 1 was rated “cloud-like” (cl), so that in terms of fit quality, this cloud is in 
about the upper one-third. 

Note that the magnetic field strength during the cloud interval is considerably 
higher than typical interplanetary values at 1 AU. Although there is no dramatic 
change in <f>, the angle 9 rotates steadily from pointing below to pointing above 
the ecliptic. The internal shock occurs slightly before 1800 UT and manifests as 
a significant jump in the magnetic field magnitude although the angles 6 and <f> 
remain largely unaffected. 

The SWE density, thermal speed, and solar wind speeds are shown in the 
bottom three panels [Ogilvie et al, 1995]. At the time of the internal shock, the 
density, thermal speed, and solar wind speed all abruptly increase, so that this 
internal shock has the characteristics of a typical fast forward MHD interplanetary 
shock. Lepping et al [1997] refer to this as a “shock-in- formation,” a compression 
still steepening into a shock, or a shock-like structure which seems to compress the 
magnetic field. We will simply refer to it as a shock. The magnetic field direction, 
as illustrated in Fig. 1, shows little change as the shock ramp passes Wind. Hence, 
it appears to be a perpendicular shock. The field is undergoing simple compression. 

In addition, inside this magnetic cloud and just downstream of the interplan- 
etary shock an unusually monochromatic wave of about one second period was 
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observed. Data from the Wind/M ASS instrument show that at the time of this 
wave there was a very cold 5 keV proton beam present, and energetic electron pitch 
angle data from the Wind/S WE instrument show that all these unusual observa- 
tions were associated with a topology change. This time period within the October 
18-19, 1995 magnetic cloud is the subject of this paper. 

A similar type of wave activity has been observed by Lucek and Balogh [1997] 
in Ulysses data. However, the period of these waves was about three seconds, and 
they lasted for about 3.5 minutes whereas the wave discussed herein had a period 
of about one second and lasted for about thirty seconds. The resolution of the 
Ulysses plasma data, however, was four minutes [Bame et a/., 1992], and although 
high resolution particle measurements are essential for a detailed analysis of shocks 
and other structures, plasma data could not be included in their analysis. 

3, Overview of Paper 

This paper is organized as follows: Section 4 gives a cursory overview of re- 
connection theory. Section 5 addresses the MFI magnetic field data describing the 
shock, upstream and downstream wave activity, and its propagation direction. Sec- 
tion 6 details the results of a variance analysis applied to the monochromatic wave 
data. Here we discuss the monochromatic wave polarization and propagation di- 
rection. Section 7 introduces Wind/3DP plasma data and shows a clear Doppler 
shift associated with this monochromatic wave. The Doppler shift allows a deter- 
mination of the wave plasma frame frequency and the magnitude of the k (wave) 
vector. Section 8 describes Wind/M ASS data which show that associated with 
this shock was a cold 5 keV proton beam. Section 9 discusses anisotropy informa- 
tion obtained from the MASS instrument. Section 10 shows that the Wind/SWE 
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instrument observed an electron heat flux drop out associated with the shock in- 
dicating a topology change. Section 11 addresses some issues dealing with shock 
orientation and propagation. Section 12 interprets the observations as resulting 
from the footpoint of the magnetic cloud reconnecting close to the Sun. Section 13 
uses the Wind/M ASS observations along with a simple coronal model to place the 
reconnection site between 1 and 5 solar radii. Finally, Section 14 provides a brief 
conclusion. 

4. Reconnection Theory 

Magnetic reconnection, a process believed to operate in and around planetary 
magnetospheres as well as in the solar corona, releases stored magnetic energy 
in the form of high velocity streams of ions and electrons and also heats up the 
particles. Figure 2 shows a simplified diagram of reconnection at an x-type neutral 
line [Kivelson and Russell, 1995]. This process is believed to result from resistivity 
breaking the frozen-in constraint within a small “dissipation region” (1) which forms 
around the x-line as shown in Fig. 2. Magnetic field lines and plasma enter the 
diffusion region from the side (2) and leave from the top and bottom effecting a 
topology change (3) in which the field lines connect to different partners. Within 
the diffusion region, the ions are accelerated away from the x-line (4), reaching the 
Alfven speed in the outflow region [Shay et al., 1999]. 

In the Petschek [1964] solution, the acceleration occurs as the plasma passes 
through slow mode shock waves connected to the diffusion region. This innovation 
increased the magnetic reconnection rate to realistic levels. A further refinement 
was made by Sonnerup [1970] who introduced fast mode shocks to loosen the inflow 
speed constraint on the Petschek model [ Kivelson and Russell, 1995]. 
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Thus, for the purpose of this work, reconnection carries with it a number of sig- 
natures: (i) shocks associated with the reconnection process (ii) topological changes 
which are frequently deduced by changing particle populations associated with dif- 
ferent flux tubes [Gosling et al, 1990]. Here, topological changes are indicated by 
the presence or absence of field-aligned and anti-field-aligned halo electrons and (iii) 
ions preferentially accelerated to the Alfven speed [Kessel et al., 1996]. 

5. MFI Magnetic Field Data 

Figure 3 shows MFI magnetic field data late in hour seventeen of day 292 (Oc- 
tober 19) 1995 around the time of the internal shock. The data are high resolution; 
MFI supplies about eleven magnetic field vectors every second at this time. The 
shock feature, between hour 17.85118 and 17.85121, makes its transition in less 
than a tenth of a second. Prior to the transition, in the upstream region, there is 
pronounced high frequency wave activity, while after the transition, there are low 
frequency, nearly monochromatic waves which, at first sight, may appear to be an 
instrumental or analysis-related artifact. 

Figure 4 shows the magnetic field magnitude in higher resolution during the 
time period immediately following the shock transition. The amplitude of the wave 
starts out at close to one nanotesla, but decays quickly over the course of about 
thirty seconds. Some residual activity is apparent for at least another thirty seconds 
after this plot, although at significantly reduced levels. This is common for the 
Earth’s fast mode bow shock due to waves being generated by thermalizing ions in 
the shock ramp. 

Using data from Geotail which observed parts of this magnetic cloud when the 
bow shock oscillated past the spacecraft, Lepping et al. [1997] found a propaga- 
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tion direction for this internal shock which is within about 20° of the cloud axis. 
Consequently, it appears that this shock is actually traveling along the magnetic 
cloud, rather than through it, suggestive of a “ducting” mode with the origin at the 
footpoints of the cloud. 

6, Variance Analysis 

Figure 5 is a hodogram showing the maximum versus intermediate component 
of the magnetic field obtained from a variance analysis [Sonnerup and Cahill , 1967] 
performed on the “first” 25 data points (about 2.5 seconds) of the wave. Five data 
points immediately after the wave rose were eliminated because they “wandered”, 
perhaps because the wave had not completely developed by the time of these five 
data points. Fig. 5 shows the wave to be elliptically polarized, as one might in 
general expect. 

The minimum variance direction obtained from the variance analysis which we 
interpret as the propagation direction is 

n = ( 0.565 , -0.599 , -0.567 ), (3) 

where the three components are the GSE x, y, and z directions, respectively, and 
the field direction is 


b = ( 0.147, -0.671, 0.726), (4) 

so that b • n = 0.073 (8 = 85.8°) so that the wave appears to be propagating 
almost perpendicular to the magnetic field (consistent with the bottom two panels 
of Fig. 3 which show little direction change in the magnetic field), but parallel to 
the shock normal (within 11°) and parallel to the magnetic cloud axis (within 30°). 
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The eigenvalue ratios were A 3 /A 1 = 11.7 and A 2 /A 1 = 2.0 so that the propagation 
direction is not well-determined, although the inferred propagation direction fits 
into our interpretation reasonably well as this propagation direction is consistent 
with the shock being the source of this monochromatic wave. 

Table 1 summarizes the relevant structure angles based on analysis from both 
this work and Lepping et al. [1997]. To a reasonable approximation, the shock 
normal, magnetic cloud axis, and wave propagation directions are aligned and per- 
pendicular to the magnetic field direction. 

Thus, this is a possible electrostatic wave propagating almost perpendicular to 
the ambient field. Such waves may be important for filling the quasi-linear gap at 
90°; that is, they will efficiently scatter particles with large pitch angles [Kammabadi 
et al , 1992]. 


7. 3DP Plasma Data 

In this study, we use high resolution plasma data from the 3DP instrument [Lin 
et al, 1995]. Specifically, the most appropriate data set is proton velocity which 
is obtained at about three second resolution. These solar wind speeds have been 
compared to the SWE and SMS/MASS data, and, although the time resolution of 
the latter two instruments is significantly lower than 3DP, the level of agreement 
is reasonably good, to about 5%. This level of agreement in the absolute value of 
the solar wind speed provides great confidence in the relative changes in solar wind 
speed observed by 3DP and used in this study. 

Because of the monochromatic nature of this wave, Fourier decomposition was 
unnecessary. Instead, to determine wave period a simple sine wave was conveniently 
fit to the magnetic field magnitude and described the data well. Fits were done 
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during the first 23.1 seconds of magnetic field data following the shock transition 
when the wave was most evident, prior to being damped out. 

Because the intention was to compare the magnetic field wave data with plasma 
data from the 3DP instrument which yields a measurement every three seconds, 
three second magnetic field data intervals were selected for fitting to a sine function. 
The eight panels in Figure 6 show the fits to eight three-second segments of MFI 
high resolution data. Table 2 lists the fit parameters corresponding to each panel 
in Fig. 6. Here mi is the d.c. offset, m 2 is the amplitude, and m 3 is the frequency. 
Table 3 compares the fit frequencies to the 3DP velocities at each time. Note the 
great stability. 

Figure 7 shows a plot of the observed wave frequency versus solar wind speed, 
which, of course, is primarily in the GSE x-direction. The data are reasonably 
described by a linear relationship with a positive slope, indicating Doppler shifting 
of the observed frequency with solar wind speed. A fit to the straight line 

tt-’obs ” kv sw COS <9 kVaw U/0, (h) 

using a value for cos $kv 3W of 0.565 from equation (3) allows a determination of 
the wavelength A = 2 tt / k and the frequency in the solar wind rest frame uq — 
u)q/'2t:. The fit shown in Fig. 7 indicates that k = 0.19 ± 0.07 km -1 . This implies 
a wavelength A = 33 km which is relatively close to the shock ramp thickness 
determined by Lepping et al. [1997] of 62 km. The plasma frame frequency u,’ 0 = 
39.7 ± 17.7 rads/s (i/ 0 = 6.3 s -1 ). As is apparent from an inspection of Fig. 7, 
there is significant uncertainty on both of these values. These values imply a phase 
velocity of u.’o/k = 209 km/s. 
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8. MASS Particle Beam Observations 

The high resolution electrostatic mass spectrometer MASS, part of the SMS 
(SWICS/MASS/STICS) package on the Wind spacecraft, was designed so that 
the particles pass through a spherical deflection system prior to entering the time- 
of-flight assembly [ Gloeckler , 1990; Gloeckler et aL , 1995; Hamilton et aL, 1990]. 
Consequently, by using start signal counts, the MASS instrument can function 
as a standard energy per charge analyzer and can determine solar wind speeds, 
proton and alpha particle densities and temperatures, and superthermal particle 
characteristics [Collier et aL , 1996; Collier et aL, 1998]. Over each spacecraft 
spin period (about three seconds) the voltage on the deflection plates changes to 
scan sixty values, so that an entire scan takes about three minutes and covers 
logarithmically an energy per charge range of 0.52 through 9.89 keV/e with a 4% 
passband. 

During the approximately thirty seconds when the wave intensity was greatest 
(1751:06-1751:33), the MASS instrument scanned the energy per charge range from 
about 3.5 to 5.4 keV/e. Figure 8 shows the background adjusted FSR2 (Front 
Secondary Electron Detection Assembly Rate-2) count rate which is unsectored, 
that is, contains no directional information. 

A statistically significant peak occurs at an energy per charge of about 5 keV/e. 
If it is assumed that this peak represents a minor ion convecting at the same speed 
as helium (minor ions tend to convect at the helium speed rather than the proton 
speed), then the species mass per charge would be about 5.83. This could correspond 
roughly to iron with charge state +10 (M/Q=56/10=5.6), not an atypical charge 
state for iron in the solar wind [e.g., Gloeckler et aL, 1999]. Although elements 
and, in particular, iron with unusual charge state distributions have been observed 
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in magnetic clouds [Gloeckler, 1998; Burlaga et aL, 1998; Skoug et aL, 1999], it 
is highly unlikely that iron of only this charge state would be enhanced and not 
the adjacent charge states, +9 and +11. Fig. 8 shows their conspicuous absence 
(M/Q=56/9=6.22 and 56/11=5.09, as shown on the upper x-axis). 

Consequently, we propose that this enhancement is most likely a 5 keV proton 
beam accelerated by a reconnection process close to the Sun. It is well-known that 
a two species, three component plasma composed of thermal ions and electrons and 
an ion beam, is unstable and will lead to wave growth [e.g. Gary et aL, 1984]. 
It is also well-known that upstream of shocks “backstreaming” and “diffuse” ion 
populations are observed [e.g. Ipavich et aL, 1984; Mobius et aL, 1986]. 

The two-paneled Figure 9 shows, a 5 a function of time, the unsectored FSR2 
count rate for the 4.91 keV/e bin (top) and the 5.16 keV/e bin (bottom). The 
highest count rate the 4.91 keV/e channel attained over the entire 28 hour magnetic 
cloud was during the thirty second monochromatic wave period. Note that (i) the 
statistical significance of the beam enhancement in the top figure is at the 3a-4cr 
level and (ii) the beam is very “narrow”. It has a low temperature and a high Mach 
number, that is the relative width in energy is less than the instrument passband 
width (AE/E < 0.051). This determines the lower limit on the velocity over the 
velocity spread, or the thermal Mach number (v/Av > 2/0.051 = 39.3). 

Finally, these ion observations contain only the energy range scanned by the 
MASS instrument over the thirty or so seconds during which the wave occurred 
(spanning 3.46-5.43 keV/e). This is a small fraction of the total energy range MASS 
scans every three or so minutes (0.52-9.89 keV/e). Consequently, the possibility of 
other ion beams or interesting features outside the observed 3.46-5.43 keV/e band 
cannot be ruled out. 
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Lutsenko and Kudela [1999] have reported more than 200 cases of “almost 
monoenergetic ions" of very short duration, of the order of one minute, with en- 
ergies between 30 and 600 keV associated with spacecraft connection to the bow 
shock. Although we do not believe the monoenergetic beam we observe has a sim- 
ilar physical origin (because the beam we observe is narrower and lower in energy, 
the wave is only on one side of the shock, and there was no change in the magnetic 
field direction across the shock), their observations are important because they hint 
that, contrary to the predictions of the standard shock associated particle models 
[e.g., Ipavich et al . , 1981 and Lee } 1982], shocks may produce high energy nearly 
monochromatic beams. 

One possible interpretation of the data involves an ~5 kV shock potential. The 
higher frequency (upstream) “hashy” waves could then be created by reflected ions 
off the potential. These reflected ions would in turn create a counter, unstable 
stream in the solar wind. Also, the energy of this counterstream would be lower 
(<5 keV), thereby resonating with the higher frequency waves and explaining why 
the emission appears at higher frequency. The low temperature of the beam may 
explain the monochromatic nature of the wave. However, the shock potential is 
generally some fraction of the upstream ion ram kinetic energy [Mandt and Kan ? 
1991], so that 5 kV may be excessively large. Thus, reconnection appears to be the 
most reasonable explanation for this monoenergetic beam. 

9. Monoenergetic Beam Anisotropies 

Figure 10 shows, during an hour interval which includes the monochromatic 
wave period, the sectored Front Secondary Electron Detection Assembly rate for 
the 45 degree sun sector, which observes flow approximately along the negative 
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GSE x-clirection, and for the 315 degree non-sun sector, which is sensitive to parti- 
cles coming from all other directions. Not surprising, there is a slight background 
anisotropy (ai ~ 0.089 ± 0.062) in approximately the solar wind direction, which 
may be simply a Compton-Getting effect due to the solar wind flow. However, 
at the time of the monochromatic wave there is a dramatic increase in the parti- 
cle anisotropy (ai = 0.363), indicating a strong anisotropy in the 45° sun sector, 
centered on the Sun-Earth line. 

Particles streaming parallel to the shock normal, which is oriented at about 
128° with respect to the GSE x-axis, would all be observed at about 52° off the 
Sun-Earth line and hence in the 315° MASS non-sun sector. Because these particles 
are observed primarily in the Sun sector and the magnetic field is primarily in the 
GSE y-z plane, we conclude that these particles have relatively large pitch angles. 
However, these anisotropy results are based on only two sectors whose averages 
could alias fine structure in the angular distributions. 

10. SWE Electron Heat Flux Measurements 

The presence of energetic electrons greater than about 100 eV or so streaming 
parallel and/or antiparallel to the magnetic field is generally viewed as a signature 
of footpoints connected to the Sun, the putative source of the electron heat flux 
[Larson et ai , 1997]. During the course of the magnetic cloud, WIND experienced 
a variety of different topologies from closed on both ends, to closed on one end, to 
open on both ends with each open field line region presumably associated with a 
reconnection event [Janoo et al. , 1998; Crooker et al, 1998; Gosling et al., 1995]. 

Larson et ai [1997] interpret the 3DP electron data during this cloud passage 
as evidence for patchy disconnection of one or both ends of the cloud magnetic field 
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lines from the Sun. The analysis of Larson et al , [1997] is not alone in suggesting 
reconnection inside the cloud. Janoo et al 's [1998] results may also be interpreted 
in this manner. 

Figure 11 shows WIND/S WE electron pitch angle distributions during a twenty 
minute period including the time of the monochromatic wave at four energies, 94 eV, 
139 eV, 203 eV, and 298 eV. Prior to the time of the shock, the electrons appear 
to be streaming antiparallel to the magnetic field indicating connection to the Sun 
on one end of the magnetic cloud. The shock appears to broaden the pitch angle 
distribution somewhat, and about nine minutes later, the 180° pitch angle electrons 
disappear indicating that the topology of the magnetic cloud has transitioned from 
being connected on one end to being disconnected on both ends, presumably due to 
a reconnection event. In fact, this may be the strongest argument for reconnection 
associated with this internal shock. As further evidence supporting this interpreta- 
tion, the internal shock is moving away from the recently reconnected side of the 
cloud, as expected if reconnection is the source of the disturbance. 

Most observed changes in cloud topology evidenced by the electron pitch angle 
distributions are interpreted as due to the spacecraft moving between flux tubes 
with different topologies and are not associated with the magnetic field, plasma, 
and energetic particle signatures seen around 1750 in the 19 October 1995 cloud. In 
order to observe these reconnection signatures, the spacecraft must be on the correct 
field line at whatever time the shock passes. Before the time near shock passage, 
the field lines will appear connected, and following the time near shock passage, the 
field lines will appear disconnected. There is a short “window of opportunity,” but 
if this interpretation is correct, eventually other magnetic clouds should be found 
which show similar shock or shock-like structures. 
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In fact, a cursory examination of magnetic field data from 34 magnetic clouds 
observed by Wind between 1995 and 199S has found three cases, including the 
October 1995 case discussed here, of shocks internal to magnetic clouds. So, it 
appears that internal shocks may manifest in about 10% or so of magnetic clouds 
observed at 1 AU. 

In addition, the magnetometer on Ulysses appears to have observed a shock 
internal to a magnetic cloud at about 5 AU late on day 228 in 1997 [Forsyth et al, 
1999] so that this phenomenon is not restricted to 1 AU. 

11, Shock Orientation 

One of the unusual features of this internal shock is its orientation which is 
roughly perpendicular to the magnetic field and anti-parallel to the cloud axis (see 
Table 1). Chao et al. [1999] have suggested that the origin of this internal shock is 
an x-ray flare located at N09W54 on October 16, 1995 at 1221 UT and claim that 
the interaction of this solar interplanetary disturbance with the magnetic cloud 
would produce the observed shock orientation. 

However, for a number of reasons, we believe reconnection near the footpoints 
of the CME to be a more natural explanation for the observations as a whole as 
well as for the shock orientation: (1) The electron observations indicate a topology 
change which is not predicted by the mechanism proposed by Chao et al. (2) The 
direction of the shock motion is consistent with an origin at the side of the cloud 
attached to the Sun initially, as deduced from the electron observations (although 
the direction of shock motion is also consistent with the Chao et al. scenario). (3) 
The energetic ion observations are consistent with reconnection close to the solar 
surface, which is not predicted by the mechanism proposed by Chao et al. (4) 
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The GOES data for the flare which Chao et al. propose as effecting the internal 
shock is extremely small, barely above background, and may not actually qualify 
by post-SOHO standards as a “flare”. Also, it could be argued that flares are not 
the best way to look for CMEs which would drive a shock, and that frequently 
very small flares can have CMEs with them. So, as a further check, SXT data 
w r ere examined [ Nariaki Nitta, private communication], but no evidence suggesting 
a CME was found. Finally, we consulted the Mauna Loa data [Joan Burkepile, 
private communication] for any reports of CMEs. Their white light observations on 
16 October 1995 began at about 1726, five hours following the flare at 1221, and it 
is quite possible that five hours after the flare a CME would be visible in the Mauna 
Loa field of view. Nevertheless, their observing logs report u no coronal activity’ 1 
for that day. In summary, although a small flare occurred on 16 October 1995, we 
have been unable to find any conclusive evidence whatsoever for a CME, let alone 
a CME driving a shock, which is necessary for the Chao et al. mechanism. 

Finally, because the main motivation of the proposal of Chao et al., was to 
explain this internal shock’s unusual shock normal direction (56° with the negative 
GSE x-axis) by postulating an additional interacting shock off the side of the cloud, 
it should be pointed out that because we expect magnetic clouds to “duct” fast-mode 
waves, the shock orientation is also consistent with the reconnection interpretation, 
as stated earlier. 

Wave ducting is a well-accepted process for coronal loops close to the Sun 
[Ireland, 1996; Smith et ah, 1997]. In the magnetic cloud, the fast mode speed 
is significantly higher inside the cloud than outside due to its higher magnetic 
field strength and lower density, so the phase fronts will move quickly within the 
cloud. Furthermore, the flux rope geometry provides that propagation near the 
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boundary is frequently quasi-perpendicular to the magnetic field and hence travels 
faster than the parallel direction [Hu, 1998]. Finally, dissipative processes such 
as viscous dissipation which prevent distrubances from propagating far from their 
source may be most effective in weak field regions [Roberts et al, , 1984; Gordon and 
Hollweg, 1983]. It may be the case that only disturbances in strong field regions 
(such as in a magnetic cloud) propagate far from their point of origin [ McLean et 
al, 1971]. 

One can determine whether or not a locomotive is approaching by placing an 
ear to the railroad tracks. Because sounds travels faster in the metal tracks, the 
train can be heard in the tracks but not in the air. It is our contention that the 
magnetic cloud plays the role of the railroad tracks, allowing us to “hear” the remote 
reconnection event when we would not be able to detect the signatures outside the 
cloud. 

12, Reconnection Model 

The reconnection interpretation of these observations is summarized by Fig- 
ure 12 which illustrates a few field lines (red, black and blue) in the magnetic cloud. 
According to the electron pitch angle data, Fig. 11, which shows anti-field aligned 
energetic electron fluxes coming only from one direction, initially one leg of the field 
is attached to the solar surface while the other has already reconnected. Then, (1) 
field lines close to the solar surface reconnect and are flung outwards at the local 
Alfven speed, about 1000 km/s (the speed of the proton beam observed by MASS), 
along with the particles on the field lines at the time. This explains the MASS 
observations of ~5 keV protons associated with the shock-they reflect the Alfven 
speed at the reconnection site. This beam may also supply the free energy for the 
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monochromatic wave. The reconnected field lines initially move much more rapidly 
than the solar wind, accelerating the solar wind downstream of the disturbance and 
creating the slight increase in solar wind speed associated with the increase in mag- 
netic field strength. (2) Following the reconnection process, the halo electrons can 
no longer populate the field lines. (3) The reconnection creates a disturbance in the 
field lines which propagates mainly within the cloud (because of the reasons cited 
in Section 11) at the local fast mode speed (the shock velocity in the solar wind 
frame is 137 km/s and the Alfven speed is 120 km/s). Thus, the shock observed 
within the cloud is propagating approximately along the cloud axis (see Table 1). 
(4) Because the reconnected field lines are thrust outwards at about 1000 km/s, and 
the field lines in front of them are travelling more slowly, the field after the shock 
passes maintains the same direction, but is compressed and increases in magnitude 
as observed in the data shown in Fig. 3. (5) The lack of field-aligned halo electrons 
suggests that throughout this process the other end of the magnetic cloud was not 
attached to the Sun. 


13. Coronal Model and Reconnection Location 


If this reconnection interpretation is correct, then the observed proton beam 
speed of about 1000 km/s represents approximately the Alfven speed at the point 
of reconnection. To determine the range of distances from the solar surface such 
an Alfven speed may mark, we examine a simple isothermal hydrostatic two fluid 
coronal model described by 


T 


dp 

dr 


GM©m p 


p + eE, 


( 6 ) 


for the protons and 



= — eE, 


(7) 
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for the electrons. Here, T is the temperature, p is the density, G is the universal 
gravitational constant, M 0 is the solar mass, m p is the proton mass, and r is the 
distance from the Sun. Although taking the electron and proton temperatures 
to be identical constitutes a bad assumption in the solar wind, for the purposes of 
determining the Alfven speed it is not critical because the mass density is determined 
by the protons and the scale height by the average of the proton and electron 
temperatures. 

Subtracting the two equations allows a determination of the electric field 


E = 


GM 0 m p 
2 r 2 e 


P- 


( 8 ) 


Plugging this back into the ion equation (4) and solving for the density yields 

. . f 1.15xl0 7 , r 0 x 1 

p(r) = ^ 0 exp| — — (l- — )j, (9) 

where T|k is the temperature in Kelvin. If the expansion is assumed radial and flux 
is conserved, then the magnetic field magnitude may be determined by 


B(r) = Bo( — ) 2 , 

r 


( 10 ) 


and the Alfven speed, v A , using the relation 


v A 


B 



(ii) 


where B is in Tesla, p is in kg/m 3 and v A is in m/s. 

Using this model with a reasonable range of values for T, po and B 0 at one solar 
radii produces an Alfven speed of about 1000 km/s (the speed of the proton beam 
observed by MASS) at a distance of between 1 and 5 or so solar radii. Figure 13 
shows the model results for Alfven speed and density for T = l.GxlO 6 K, po = 



Collier et al.: Reconnection Remnants in a Magnetic Cloud 21 

l.TxlO -12 kg/m 3 , and B 0 = 1.3xl0 -3 T. The 1.3xl0 -3 T magnetic field value 
represents a radial extrapolation of the observed 28 nT Wind field back to the solar 
surface. As the field probably expands superradially, this likely underestimates the 
field somewhat. Also, because the slow (i.e. equatorial) solar wind flow is believed 
to be subsonic to about four solar radii, the inclusion of solar wind flow will not 
affect the model Alfven speed significantly in the region of interest. These radial 
distances for reconnection are reasonable and support the scenario that W IND is 
observing reconnection remnants inside the October 18-19, 1995 magnetic cloud. 

14. Conclusion 

We have presented WIND observations from MFI, 3DP, SMS, and S\\ E inside 
the October 1S-20, 1995 magnetic cloud which may be interpreted as evidence of 
reconnection occurring at low altitudes in the solar corona, between 1 and 5 solar 
radii. The observations include an internal shock travelling approximately along the 
axis of the magnetic cloud, simple compression of the magnetic field consistent with 
the footprint magnetic fields being thrust outwards at speeds much greater than 
the solar wind speed, an electron heat flux drop out occurring within minutes of the 
shock indicating a topological change resulting from disconnection from the solar 
surface and a very cold 5 keV proton beam resulting from reconnection. In addition, 
an unusually monochromatic wave propagating perpendicular to the magnetic field 
was observed in association with the shock and may be related to the monoenergetic 
particle beam. 

Because the spacecraft must be fortuitously positioned to observe the shock 
when it passes, observations such as the ones reported here will be rare. Further- 
more, because shocks will expand, it may be possible to observe a shock internal 
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to a magnetic cloud without a corresponding topology change or energetic particle 
beam. However, if our interpretation is correct, given observations of enough mag- 
netic clouds with sophisticated instrumentation, similar observations inside mag- 
netic clouds should surface in the future, particularly as we transition to solar 
maximum and magnetic clouds become more frequent. Consistent with this, a few 
magnetic clouds have been found to possess internal shocks. 
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Figure Captions 

Figure 1 . Overview of the Wind MFI and SWE magnetic field and plasma obser- 
vations during the October 18-19, 1995 magnetic cloud. The interior of the cloud 
is marked with the two-headed arrow, although there is some ambiguity about the 
back of the cloud. The internal shock is most evident in the magnetic field mag- 
nitude, in the first panel, whereas the magnetic field direction, given in the second 
and third panels, varies little across the shock. The density, thermal speed and ve- 
locity all increase across this shock so that it shows characteristics of a fast forward 
MHD shock. The shock speed is very close to the Alfven speed, approximately the 
fast-mode speed. 

Figure 2. Reconnection process at an x-type neutral line. Magnetic field and 
plasma are accelerated in the diffusion region (1) after entering from the left and 
right (2). The field lines and plasma are ejected towards the top and bottom 
effecting a topology change (3). The ions are preferentially accelerated to the Alfven 
speed (4). 

Figure 3. Overview of the MFI high resolution (eleven vectors every second) 
magnetic field data showing the “shock” inside the October 18-20 magnetic cloud. 
The top panel shows the magnetic field magnitude in nanotesla with the second 
through fourth panels showing the individual GSE components of the magnetic field. 
The fifth and sixth panels show the 9 (out of the ecliptic) and (f) (in the ecliptic) 
angles of the magnetic field. The shock transition occurs between hour 17.851 IS 
and 17.85121 on day 292 ( October 19) with the monochromatic wave immediately 
following and lasting about thirty seconds. There are hashy high frequency waves 
(3>1 Hz) present on the upstream side. Note that the magnetic field angle changes 
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very little during this time period. 

Figure 4. The magnetic field magnitude during about thirty seconds immedi- 
ately following the shock transition. The monochromatic wave activity is most 
pronounced but decaying in amplitude during this thirty second interval, although 
some trace of the activity, at a much lower level, is present for a minute or so after 
the shock transition. 

Figure 5, Hodogram showing maximum versus intermediate component for the 
“first” 25 data points of the wave. The eigenvalue ratios of the variance analy- 
sis are Aj nt /A m i n = 2.0 and A max /A; nt = 5.7 so that the propagation direction is 
not exceptionally well-determined. Note that the wave appears to be elliptically 
polarized. 

Figure 6. These eight panels show the sine fits to three second intervals of MFI high 
resolution magnetic field magnitude data. The shock jump occurs at hour 17.85121 
with a wave significantly damped by about hour 17.8593. The wave duration was 
23.1 seconds and the functional form fit was mi+ni2*cos(3600*ni3 *(mo —start time)) 
where “start time'" is the beginning of the three second data interval. Phase was 
eliminated by starting and ending on a complete half cycle. The values for mj , m 2 , 
and m 3 , as well as the start time used for each panel (a-h) are listed in Table 2. 

Figure 7. The observed frequencies in the Wind spacecraft frame obtained from 
the fits shown in Fig. 6 plotted versus the observed 3DP three-second solar wind 
speed data. If the Doppler shifting is assumed to obey cj 0 bs = kv sw cos #kv 3W + 
where k is the wave vector and uj 0 is the rest frame wave frequency, then the fit given 
by the solid line provides u>q — —39.7 ± 17.7 radians/s and k = 0.19 =t 0.07 km -1 . 
The correlation coefficient for these data is 0.73. 
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Figure 8, The background adjusted SMS/MASS FSR2 count rate during the wave 
time period (17 5106-5133). The time periods used for background subtraction were 
17 4500-4528 and 17 5409-5437, which cover the same energy per charge range. A 
clear, statistically significant peak occurs at an energy per charge of about five 
keV/e. The mass per charge (indicated at the top of the plot) corresponding to 
iron of charge state +10 is indicated by the arrow, although it is argued the beam 
is probably suprathermal protons. 

Figure 9. Counts in the 4.91 keV/e and 5.16 keV/e energy per charge bins as 
a function of time over a one hour interval from 17:20-18:20. The beam is very 
“narrow” in velocity space. 

Figure 10. The sectored FSR rate for both the 45° sun and the 135° non-sun 
sectors for the energy per charge bin 4.91 keV/e are plotted on the left and right y- 
axes, respectively. The logarithmic y-axes are offset by a factor of seven to account 
for the disparate angular ranges of the two sectors. 

Figure 11. WIXD/SWE electron pitch angle distributions around the time of 
the internal shock-like feature. The four panels show four energies: 94 eV, 139 eV, 
203 eV, and 298 eV. Prior to the shock, the electrons appear to be primarily stream- 
ing at 180° pitch angle. The shock appears to broaden the pitch angle distribution, 
and about nine minutes after the shock, the streaming electrons disappear. 

Figure 12. Schematic summary of the reconnection interpretation for the WIND 
observations inside the magnetic cloud of October 18-19, 1995. Initially, a field 
line disconnects from the solar surface and is flung outwards at ~-T000 km/s. The 
reconnection event triggers a disturbance in the field lines which propagates along 
the cloud axis. The field lines before and after the discontinuity maintain their 
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direction, but are compressed. 

Figure 13. S ample results for the Alfven speed and the density as a function of 
distance from the solar surface for a simple isothermal hydrostatic two fluid coronal 
model using a coronal temperature of l.GxlO 6 Kelvin, a density of 1.7xl0 -12 kg/m 3 , 
and a magnetic field of 1.3xl0” 3 Tesla. This model is used to place the location of 
the reconnection point between 1 and 5 solar radii. The one solar radius field value 
used represents a radial extrapolation of the 28 nT field observed at Wind back to 
the Sun and as such represents a lower bound for the field at the Sun since most 
probably the field expands superradially. 
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Matrix of Angles 



cloud axis 

shock normal 

wave propagation 

shock normal 

159° (||) 

— 

— 

wave propagation 

31° (II) 

169" (||) 

— 

magnetic field 

59° (_L) 

105° (±) 

86° (_L) 


Table 1 
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panel 

nil (nT) 

m 2 (nT) 

m 3 (rads/s) 

start time (hr) 

a 

27.99 

0.66 

6.03 

17.85218 

b 

27. SO 

0.75 

6.10 

17.8530S 

c 

27.73 

-0.63 

5.98 

17.853S4 

a 

27.71 

-0.49 

6.37 

17.85471 

e 

27.68 

-0.62 

6.46 

17.85558 

f 

27.65 

0.49 

6.24 

17.85647 

g 

27.62 

0.46 

6.41 

17.85732 

h 

27.65 

0.25 

6.61 

17.S5816 


Table 2 
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time (hour) 

MFI ljq (rads/s) 

3DP | v | (km/s) 

17.85256 

6.03 ± 0.03 

424.28 

17.85341 

6.10 ± 0.05 

421.83 

17.85426 

5.98 ±0.05 

424.02 

17.85510 

6.37 ±0.03 

425.93 

17.85595 

6.46 ±0.04 

426.32 

17.85680 

6.24 ± 0.04 

424.31 

17.85765 

6.41 ±0.05 

424.31 

17.85850 

6.61 ±0.06 

426.32 


Table 3 
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